Astronomy & Astrophysics manuscript no. main 


©ESO 2011 


October 5, 2011 
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ABSTRACT 

Context. Recent observations of stars with unusually large radial velocities in the Galactic stellar halo have raised new interest on 
so-called Hypervelocity stars. Traditionally, it is assumed that the velocities of these stars could only originate from an interaction 
with the supermassive black hole in the Galactic center. It was suggested that stars stripped-off a disrupted satellite galaxy could reach 
similar velocities and leave the Galaxy. 

Aims. In this work we study in detail the kinematics of tidal debris stars to investigate the implications of the new scenario and the 
probability that the observed sample of Hypervelocity stars could partly originate from such a galaxy collision. 
Methods. We use a suite of Af-body simulations following the encounter of a satellite galaxy with its Milky Way-type host galaxy to 
gather statistics on the properties of stripped-off stars. We study especially the orbital energy distribution of this population. 
Results. We quantify the typical pattern in angular and phase space formed by the debris stars. We further develop a simple stripping 
model predicting the kinematics of stripped-off stars. We show that the distribution of orbital energies in the tidal debris has a 
typical form which can be described quite accurately by a simple function. Based on this we develop a method to predict the energy 
distribution which allows us to evaluate the significance and the implications of high velocity stars in satellite tidal debris. 
Conclusions. Generally tidal collisions of satellite galaxy produce stars which escape into the intragalactic space even if the satellite 
itself is on a bound orbit. The main parameters determining the maximum energy kick a tidal debris star can get is the initial mass of 
the satellite and only to a lower extent its orbit. Main contributors to an unbound stellar population created in this way are massive 
satellites (M sat > 10 9 M Q ). We thus expect intragalactic stars to have a metallicity higher than the surviving satellite population of 
the Milky Way. However, the probability that the observed HVS population is significantly contaminated by tidal debris stars appears 
small in the light of our results. 
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1. Introduction 

The growth of galaxies via the accretion of smaller companion 
systems is one of the major ingredients in the current percep- 
tion of galaxy formation and evolution. These satellite galaxies 
are disrupted in the tidal field of their host galaxies and the new 
material is dispersed near the orbit of the progenitor. Recent the- 
oretical work has shown that especially the outer stellar halo is 
predominantly made of stars which were born outside the main 



galaxy ( Abadi et al. 2006 ; Zolotov et al. 2009 ; Scanna pieco et al.| 
|201 l| i. Such stars are thought to have low metalhcities and to be 
old. The small fraction of material born inside the main galaxy 
reaching these large radii was mostly re-distributed during vio- 
lent major merger events. As these events were more frequent 
when the galaxy was still young, these stars are also predomi- 
nantly old. A third small population of the outer halo are the so- 
called Hypervelocity stars (HVSs) which are ejected via a three- 
body interaction involving the supermassive black hole (SMBH) 
in the Galactic center ( |Hills|19 88). Such stars have no age con- 
straints and should be metal-rich as they originate from the in- 
nermost region of the galaxy. 

The latter population earned attention since they could serve 



as an indirect proof for the SMBH in the Galactic center (Hills 
1988) and also because they could be used to measure the shape 
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of the Galactic po t ential (|Gnedin et al.|2005[|Yu"& Madau|2007 



|Perets et al.|2009) . |Yu & tremaine|(|2063D estim ated a HVS ejec 
tion rate of 10~ J yr -1 and |Perets et aT 72007 1 showed that this 



rate could increase by a magnitude if massive perturbers such as 
giant molecular clouds or star clusters were considered. Aside 



from the classical ejection mechanism by Hills ( 1988 1 several 



authors have suggested alternative formation scenarios: a binary 
black hole of equal ( Yu & Tremai ne|2003 1 and un-equal masses 
(Levin 2006; Sesana et al.||2009| l or the accretion of a satellite 



galaxy ( |Abadiet al. 2009 



Recent observations of 17 stars in the Galactic halo with un- 



2005 



Hirsch lTaT1[2005 



b| |2007a|b| 2009a 



usually high velocities ( Brown et al 
Edel mann eFaLl[2003] |Brown et al.||2006a 
Tilli ch et al.||2009) > raised new interest on the topic. By de- 
sign of the search strategy these stars have typically blue colors. 
They move with velocities up to 720 km s _1 with respect to the 
Galactic center and are thought to reside at Galactocentric dis- 
tances of 20-130 kpc. Interestingly, the targeted HVS survey of 



|Brown et al. ( |2009a| l only yield out-going HVSs which is typ- 
ically attributed to the short lifetimes of the stars compared to 
the long orbital periods. However, eventually also an in-going 
star with extremely high velocity was observed (Przybilla et al. 
|20T05] >. 

Despite the small number of HVSs reported to date several pe- 
culiarities in the distribution of the observed population were al- 
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ready claimed. fAbadi et aT1 p009 ) found that a large part of the 
population clusters around a certain travel time (~ 130 Myr), i.e. 
the time a star would need to travel from the GC and arrive at its 
current radius with its current radial velocity. Such a clustering 
is not expected from Hills' original SMBH-ejection scenario. It 
could be explained, however, by a star burst event near the GC 
triggering an increased ejection rate of HVSs for certain times. 
Also the angular distribution on the sky shows signs of 
anisotropy. |Abadi et al.| ( |20Q9] ) and Brown et al. ( 2009b) > report a 
significant overdensity of HVSs around the constellation of Leo. 
Stars ejected by one or more black holes in the GC should appear 
on the sky in an approximately homogeneous distribution or in 
a ring-like structure (Levin 2006 ). However, a preferred ejection 
directory as found in the data is not naturally explained with this 
mechanism (however, see Lu et al. (2010)). 

The accreted population of stars in the outer halo can also 
contain stars with large radial velocities. Teyss ier et al.| f2009) 
showed that there should exist an energetically loosely or un- 
bound population of stars originating from disrupted dwarf 
galaxies. Aba di et aT] ( p0 09 ) commented that a larger total mass 
of the Galaxy would allow the normal virialized halo population 
to reach these velocity regimes. The authors further suggested 
that the peculiarities in the HVS distribution would be naturally 
explained if part of the observed HVSs would actually belong to 
a stream of tidal debris of a recently accreted dwarf galaxy. An 
example for a HVS likely being generated by this mechanism 
was recently found in M31 (Cald well et al.|20T0] >. 

Several theoretical studies have already investigated proper- 
ties of the tidal debris of satellite galaxies. |Johnston] ( | 1998[ > ap- 
proximated the energy distribution of tidal debris particles with 
a triangular shape to build up a stellar halo distribution. Choi 
et al.| ( |2007] > showed that the energy kick obtained by stripped 
stars via tidal forces and also the deviations between leading and 
trailing tidal arms are both increasing with the mass of the ap- 
proaching satellite. Warnick et al. ( 2008 1 investigated the rela- 
tion between observable properties of tidal streams like radial 
velocity dispersion or thickness to the properties of the progen- 
itor system. D'Ong hia et al.| ( f2"009||20l6") l investigated the effect 
of resonances during tidal stripping of rotating systems. 

In the present work we investigate the kinematical proper- 
ties of tidal debris with a special focus on the fastest stars of this 
population. For this we systematically study tidal encounters of 
satellite galaxies with their hosts and comprehensively investi- 
gate the process. We ran a suite of collisionless A^-body simula- 
tions following the passage of a small companion galaxy through 
its massive Milky Way-like host galaxy. The set-up of these sim- 
ulations is described in Sect. [2] In Sect. [3] we analyze these sim- 
ulations to obtain an idea of what properties an observer would 
find in an HVS population generated in a tidal collision. Then a 
simple analytical model is developed and tested against the sim- 
ulations (Sect. [4j. We show in Sect. [5] how this model can be 
used to predict the energy distribution of the tidal debris star. In 
Sect. [6]we discuss our results and summarize them in Sect|7] 



2. Simulation set-up 

All simulations were run using the publicly available simula- 
tion code Gadget-2 (Springel 2003} . For the main galaxy we 
used the model parameters proposed by |Klypin et al.| (|2002 ). The 
galaxy consists of three components, an adiabatically contracted 
spherical NFW halo ( Navarro et al.|1997) , an exponential stellar 
disk and a spherical stellar bulge with a Hernquist density profile 




Virial mass M M [10 10 M s ] 

Fig. 1. Black symbols show the applied values for the halo con- 
centration parameter c for a given satellite mass M sat . The er- 
rorbars have a width of 0.5 to reflect the fact that only integer 
numbers were used for c. The red line shows a power law with 
exponent 0.36. 



( Hernquist|1990 i. The disk density profile is 



Pdisk(fl,z) 



M disk 



exp 



R 

Rdisk 



sech 



\Zdisky 



(1) 



where R = y* 2 + y 2 . The vertical scale height is set as Zdisk = 
0.2Rdisk, as was found in observations of other disk galaxies 
(Kre gel et al.|2()02). Table [T]presents the parameters used. 

The satellite galaxies are modeled as an adiabatically con- 
tracted NFW halo hosting a Hernquist sphere as the spherical 
baryonic component. The models fulfill the constraints from the 
Fundamental plane of dE+dSph galaxies given in de Rijcke et al. 
( |2005) : 



(2) 
(3) 



log Lg ~ 4.39 + 2.55 log cr 
log Lg ~ 8.69 + 3.55 log R e 

Here, Lg is the fi-band luminosity, R e is the effective radius en- 
closing half the light of the galaxy and cr is the luminosity- 
weighted mean velocity dispersion. A Hernquist sphere has an 
effective radius R e =s 1.82rb u i ge (| Hernquist]|1990[ ). For the ve- 
locity dispersion <xo, in analogy to the dispersion used by de 
Rijc ke et al.| ( 2005 ), we computed the mass-weighted mean of 
the line-of-sight velocity dispersions in different radial annuli of 
the visible component. Finally, to relate the luminosity Lb to the 
baryonic mass content of the galaxy we assume a mass-to-light 
ratio T B = 2T B , . 

The ratio between total, M sat , and baryonic mass, M satj b, was 
fixed using the relation found by McGaugh et al. ( 20 1 0) > : 



logM sat , b = 4.01ogv circ + 1.65. 



(4) 



To relate the circular velocity of the dark halo, v c i rc , to the mass 
of the satellite, M sat , we used the same relation as in McGaugh 
[eTaLl ( |20T0l i: 



M sat = (1.5 x 10 5 km" 3 s 3 M o )v c , i 



(5) 
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With these constraints the dwarf galaxy is fully determined by 
only one parameter. In this work we used the total mass M sat as 
a free parameter which was varied for different simulation runs. 
The requirement to match all the constraints given above fixes 
also the concentration parameter c of the satellite dark halo. The 
obtained mass-concentration relation is given in Figure [TJ It is 
best described by a power law: 



Table 1. Parameters of the host galaxy 



C a 19.3 



10 10 M P 



0.36 



(6) 



Note, that this relation was obtained by fitting our basic satel- 
lite model to the observational constraints. Our concentration 
parameter should thus not be interpreted in the original sense 
of an evolutionary sequence in the frame-work of ACDM cos- 
mology ( Navarro et al.|1997 1. 

For the initialization of the phase-space positions of the particle 
samples we followed a method outlined by Springel & White 
( 1999) which is a modified version of the method of |Hern quist 
( 1993 i which assumes Gaussian velocity distributions. The lat- 



ter leads to slight deviations from a perfect equilibrium config- 
uration. To account for this both host and satellite galaxies are 
allowed to relax for 1 Gyr in isolation before they are imple- 
mented into the actual simulations. We use a softening length of 
0.01 (0.2) kpc for the satellite star (dark matter) particles in all 
our simulations. 

Simulation time All simulations ran until the satellite reached 
the Apogalacticon after its first pericentric passage or crossed 
the virial radius of the host galaxy, /?20o- This was done because 
we wanted to study the properties of a population of tidal debris 
particles generated during a single stripping event. As we also 
have a focus on the stars escaping from the combined satellite 
and host system we chose to study only the first orbit. We expect 
this orbit to generate the largest spread in velocities as the ini- 
tial unperturbed satellite population covers the complete possible 
phase space regions. At later orbits the satellite will have lost its 
most energetic population ( |Choi et aT]2 009). Furthermore, con- 
sidering only the first orbit allows an easier comparison between 
the different simulation runs as one has the full control over the 
satellite configuration at the beginning of the orbit. 

A suite of simulations To create a suite of comparable simu- 
lations we then ran this scenario with varying initial satellite 
masses M sat = 0.1 - 2 x 10 M G and different starting posi- 
tions in the satellite phase space. In the majority of cases the 
satellite is on a bound polar orbit with respect to the host disk 
component. To determine the influence of an inclined orbit we 
also ran a couple of simulations with 0° (planar) and 45° in- 
clination angle. We found that the differences in the results for 
varying inclinations were small compared to other uncertainties. 
We thus did not consider inclination as a major parameter and 
neglect it completely. The initial angular momenta, L sat , of the 
satellites range between and 15 x 10 3 kpc km s _1 which corre- 
sponds to pericenter distances R pe „ from to 50 kpc. Table D.l 



lists the initial conditions and some of the analysis results for all 
runs. 

3. Observable properties from the simulations 

Of the 41 simulations performed for this study 23 yielded par- 
ticles with velocities higher than the local escape speed of the 
host galaxy. These particles are gravitationally unbound and are 
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the simulated equivalent to HVSs. However, in the real Milky 
Way the escape speed is still uncertain to a c onsiderable degree 



potential (Law et al 



Smith et al. 



2007; 



as neither the total mass (1 - 2xl0 12 M Q ,e.g.[ 
I Xue et al.|[2008) |Guo et al.pOlO] [Boylan-Kolchin et al.| 
|Przybilla et al.||2010b 1 nor the global shape of the gravitational 



20TT1 



2009) ) is measured accurately. Moreover, 



the asymmetry of the Galaxy introduces a direction dependency. 
Thus the value of the escape speed must not be seen as a sharp 
limiting velocity dividing bound and unbound stars. It is rather 
a characteristic value to compare to when evaluating the prob- 
ability whether a star will eventually fall back onto its host or 
not. For the dynamics of a star located within the virial radius of 
the Galaxy it makes no qualitative difference whether it is grav- 
itationally unbound. Because of this, and also to obtain better 
statistics, we will thus in this section analyze the most energetic 
0.1% of the satellite particles, i.e. the Most Energetic Particles 
(MEPs, for short) regardless whether they reach velocities higher 
than their local escape speed. 

Angular distribution Figure [2] shows the angular distribution of 
the MEPs in one of our simulations (red stars) as seen by an 
observer on the Sun's location. The snapshot is taken at a time 
when the particles have moved from the satellite Perigalacticon 
out to a galactocentric distance of ~ 60 kpc, similar to the ob- 
served HVS population. Green dots indicate the positions of all 
satellite particles while black dots represent star particles be- 
longing to the host galaxy. As already reported by Abadi et al. 
(2009) the MEPs are clustered in a tightly confined region on 
the sky (in Fig. [2] marked by the solid magenta circle of radius 
~ 15°). Averaging over all 39 simulations and over 10 equis- 
paced angular positions of the sun on the solar circle the mean 
radius of a circular region encompassing all MEP is 16° and the 
maximum angular radius is 27°. At a distance of 60 kpc the view- 
ing angle of the observer relative to the orbital plane of the satel- 
lite does not play a significant role as the stripped-off particles 
had not enough time yet to unfold into a prominent stream and 
are thus observable in a compact area from all directions (cf. the 
upper left panel of Fig. [5J. 

Furthermore, the position of the satellite relative to the MEPs is 
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-75° 



Fig. 2. Aitoff projection of a simulation run about 200 Myr after Perigalacticon. Green dots represent all satellite particles (only 
every 5th particle plotted) while the MEPs are marked by red stars. Disk and bulge of the host galaxy is also shown as black mass- 
density contour-lines. The median distance of the MEPs is ~ 60 kpc, comparable to the Galactocentric distances of the observed 
HVSs. The MEPs are concentrated to confined area on the sky, in this case enclosed by a circle with angular radius 15°. 



not arbitrary. In all our simulations the satellite has a smaller an- 
gular distance to the host galactic center then the MEPs. This 
is due to the fact that the satellite as well as the MEPs have 
just passed their perigalacticon and now move away from the 
host. Since the MEPs have higher orbital energies they leave the 
satellite behind and are thus observed at larger angular distances 
in the vast majority of cases. The angular distance between the 
MEPs and the satellite remnant is determined by projection ef- 
fects depending on the viewing angle relative to the orbital plane 
of the satellite. In our simulations the satellite core is observed 
always within 22° (90 percent within 16°) separation to the cen- 
ter of the MEP population. 



The radial velocity-distance plane The particles stripped-off 
the satellite quickly disperse in physical space and are soon 
indistinguishable from the already existing Galactic halo pop- 
ulation. However, when Galactocentric radial velocities (RVs) 
are plotted against Galactocentric distance, r, the particles form 
an elongated pattern reflecting their common origin from the 
Perigalacticon of the satellite. Figure [3] shows time series of two 
simulations with satellite systems on orbits with different ec- 
centricity. Green dots represent satellite particles while red stars 
show the MEPs. Particles of the host galaxy are not plotted. The 
upper series is based on the same simulation that was used for 
Fig. [2] and the upper left most panel shows the same point of 
time. The lower panel row shows a simulation run with an al- 
most purely radial orbit which has also lower orbital energy. The 
initial structural properties of the satellite are the same as in the 
upper run. Note that despite the lower satellite orbital energy the 
maximum tidal debris velocities are larger than in the run with a 
lower eccentricity. 

The dashed line represents the escape speed from the host sys- 
tem at the respective distance. This corresponds to the trajectory 
of a test body on a parabolic (purely radial) orbit. Particles be- 
low this line are not necessarily bound as part of their motion 
is hidden in the other velocity components. This explains why 
particles can cross the dashed line while still conserving their 
energy. 

The dotted lines in the panel are lines of constant travel time, 



i.e. mark positions in the RV-r-plane which are occupied by test 
bodies which started from a common point with varying initial 
(radial) velocities. The starting point is usually chosen to be the 
center of the gravitational potential, as in Fig. [3] However, lines 
from other starting points have similar shapes which explains 
the good alignment of the tidal debris particles even though the 
satellite in the upper row run never comes close to the Galactic 
center (the pericenter distance for this run was 18 kpc). 
The width of the tidal debris streams in the RV-r plane also 
changes with eccentricity. This can be understood when consid- 
ering that most of the stripping happens during a short period of 
time around the pericenter passage. For more circular orbit this 
period get more extended as more time is spent by the satellite 
at radii similar to the pericenter distance. 



Hence the clustering in travel time already reported by Abadi 
|et al.| ( |2009] > has an intrinsic scatter which scales with eccentric- 
ity of the orbit of the progenitor system. This is especially so 
shortly after the pericentric passage, when the distances to the 
GC are not large and lines of constant travel times have a large 
slope. At this time the scatter can be stronger than the trend to 
lines of constant travel time. 

The left more panel (upper row) shows the distribution of the ob- 
served HVSs in the distance-velocity plane for comparison. The 



possible tidal debris group in the population proposed by Abadi 
|et al.| ( |2009] l clusters around the 133 Myr-line of constant travel 
time show as dotted black line. Note, that none of the two sim- 
ulation runs was designed to reproduce the observed population 
as this was not the goal of this more general parameter study. 



Maximum velocities In the classical SMBH sling shot scenario 
the extremely large ejection velocities are a result of the extreme 
orbital velocities occurring near a SMBH plus the large orbital 
velocities of the components of a hard binary system ( |Yu &| 
Tremaine 2003). Compared to such an environment, collisions 
of galaxies are much less violent events as the time scales are 
much larger and potential gradients much shallower. We thus 
cannot expect the extraordinary velocities up to 3000 km/s pre- 
dicted by Hills ( 1988 ). Still, the simulations show that stars can 
be accelerated to their local escape speed and above. The maxi- 
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Fig. 3. Time series of the galactocentric RV-distance plane for two simulation runs. The time elapsed since Perigalacticon, Af per i, is 
shown in the upper left corner, respectively. Red stars represent the MEPs, green dots satellite particles and dashed lines the local 
escape speed of the respective radius. Dotted lines mark lines of constant travel times corresponding to Af per i. The eccentricity of the 
satellite orbit is given in the lower right corner of the right most panel, respectively. In the left most panels the mean distance of the 
MEPs is 60 kpc similar to the observed HVS population which is shown in the right most panel in the upper row for comparison. 
At this point the MEPs did not have time yet to leave the remaining satellite far behind. Note, that none of the simulations were 
designed to reproduce the observed HVS population. 



mum velocities reached by the MEPs at r a 60 kpc in our sim- 
ulations range between 200 and 400 km/s (v esc (60kpc) = 330 
km/s in our Galaxy model). 

4. A model for tidal stripping 

To guide our further analysis we develop a simple, succinct 
model to describe the accreted satellite mechanism. It was in- 



spired by the calculations of HVS ejection velocities by Yu & 
Tremaine ( 2003[ ) as it treats the galaxy-galaxy encounter similar 
to an binary-SMBH encounter.: a satellite galaxy is moving on 
an orbit in the gravitational potential <I>host( r ) of its much more 
massive host galaxy. Its specific orbital energy is thus 



+ G>host(r sat ) 



(7) 



Since the satellite is an spatially extended object it is subject to 
tidal forces which lead to a mass loss of the satellite. Under the 
assumption of an at least moderately eccentric orbit the major- 
ity of this stripping will happen in a short period of time when 
the satellite is closest to the center of the host galaxy where tidal 
torques are strongest, i.e., at its perigalacticon, /? pe ri, where it has 
the velocity v sat = V peri . 

To model the stripping we now assume what we call instanta- 
neous escape: a star i with a position r, relative to the satellite 
center and a velocity v, in the co-moving rest frame of satellite 
has an orbital energy 

1 9 

Ei = -(Vsat + v,) + <IW(r S at + r,-) + Osatfa). (8) 
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Fig. 5. Energy gain with respect to E s . dt predicted by our simpli- 
fied model AE m . dx compared to the maximum energy gain found 
in each simulation. The latter are by a factor of ~ 0.45 (slope 
of the grey dotted lines) lower than the estimates, which is most 
likely due to oversimplification of the model. Phase-space sam- 
pling due to the limited particle number does not play a signif- 
icant role as can be seen from the grey circles which mark the 
runs with 10 times lower resolution. Left: the estimated energy 
gain as obtained from our model. Right: the energy gain from 
our model corrected for an additional dependency on the angu- 
lar momentum of the satellite. See text for a discussion. 



It is lost to the satellite when the gravitational potential from 
the satellite, <D sat is less than the difference in the host potential 
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Fig. 4. Snapshots of the satellite (small green dots) shortly before, during and shortly after its pericentric passage. Red stars show the 
positions of those particles which will have the highest orbital energy at the end of the simulation. Blue triangles show those with 
the lowest energy. The latter have a relative motion opposite to the satellite motion and are on retrograde orbits. The high energy 
particles move with the satellite and are on prograde orbits. 



between the satellite position and the position of the particle, r,-, 

<t>ho S t(r sa t + r,) - O host (r sat ) > -O sat (r,), (9) 

which is equivalent for it to be outside of the tidal radius, ^ t id a i- 
We now assume that this energy transition occurs instantly and 
the star is left to move in the host potential only. Thus the orbital 
energy of the star after the stripping is 

E'i = ^(V sat + V,) 2 + $host(r S at) = 

1 , 1 , 
= V sat Vi + - vf + -v^t + O host (r sat ) = 

1 2 

= VsatVi + -V| + East. 



The energy kick obtained by the star compared to the satellite is 
then 



1 ? 

AEj = V sat V; + -v,. 



(11) 



We now ask for the maximum of AE. Equation [TT| leads to the 
assumption that three conditions need to be fulfilled for the max- 
imum energy kick: 

1 . The star has the maximum velocity possible which is the lo- 
cal escape speed at the tidal radius: v, = v e sc(^tidai)> 

2. Satellite and stellar velocities have to be aligned: v sat ||v,-, 

3. The satellite has to be at its maximum velocity, which occurs 
during the passage of the Perigalacticon: v sat = V per i. 

Moreover, if the star is to gain orbital energy it needs to be at 
larger Galactocentric radii than the satellite, i.e. |r sat + r,| > r sat , 
because only then the tidal force push the star away from the 
galactic center and thus from the potential well. This, together 
with the first condition, requires the star to be on a prograde orbit 
with respect to the satellite motion. This view is also confirmed 
by Fig. [4] It shows three snapshots of a simulation run shortly 
before, at and shortly after Perigalacticon. Red stars and blue 
triangles mark those particles which will have the highest/lowest 
orbital energy at the end of the simulation, respectively. The two 
groups are situated are very distinct locations with respect to the 



satellite. The particles which gain energy move along with the 
satellite on an orbit prograde with respect to the satellite motion 
while for particles which lose energy the orbital phase is such 
that their motion is contrary to the direction of the system veloc- 
ity. 

Thus we arrive at 



1 2 

A£m ax - Vperi V esc (^tidal) + ~ Vesc(^tidal) • 



(12) 



In the left panel of Fi g. [5] we compare this model prediction ob- 



(10) tained from Equation 



12 to simulation results. For the latter we 
analyzed our simulations at Apogalacticon after the satellite's 
passage through the host system and the particle with the highest 
orbital energy max(£',) was identified. The energy gain AE max s j m 
is then defined as max(£ , ,)-£ , sat . For further details how this esti- 
mate was done we refer the reader to the Appendix|A] Solid lines 
in Fig. [5] indicates equality. Equation 12 systematically overes- 
timates the maximum energy kick by a (constant) factor of 2.2, 
as indicated by the dashed line with slope a = 2.2~ l 0.45. 
The overestimation are most likely due to oversimplification of 
our estimate, e.g. the neglection of the structural evolution of the 
satellite in the tidal field. Surprisingly, the mass resolution only 
plays a minor role in the sense that simulation runs with lower 
particle numbers do not yield significantly smaller maximum en- 
ergy differences. This is most likely due to the steep slope of 
the energy distribution of the tidal debris stars (see Sect. [5] and 



Fig. B.l 



Much of the scatter in the left panel of Fig. [5]turned out to be a 
residual dependency on the initial angular momentum of the re- 
spective satellites in the simulation. We obtain a tighter relation 
(right panel) when we use 



L S3f + L 



char 



(13) 



with L c h al = 6300 kpc km s _1 . Thus the value aAE maxc gives a 
robust estimate of the maximum energy gain occurring during a 
satellite-host galaxy encounter. 
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Dynamical friction In the course of its orbit the satellite galaxy 
is also subject to dynamical friction. This means that it will sink 
deeper into the potential well of the host system. By the time 
it reaches its Perigalacticon stripped-ofF stars might have to get 
over an energy gap much wider to become HVSs. The extent 
of the energy loss depends mainly on the mass M sat and the or- 
bit of the satellite, in a way which is counterproductive to the 
ejection energies AE max : a more massive satellite ejects stars at 
higher energies as its internal velocities are larger. This energy 
gain for the stripped-off stars goes roughly spoken with y/M*,,, 
since v esc (r) oc y/M sat which goes into our estimate (Eq. 12 1. On 



the other hand a higher mass results in stronger dynamical fric- 
tion which roughly goes oc M sat (Chandrasekhar 1943). In fact in 
our simulations we find that the loss in orbital energy after one 
orbit is 



A£ DF = 2 x 10" 



Msat 



0.7S 



1 + 



5900 kpc kms" 



(14) 



which reflects the dependency of the physical extent of the satel- 
lite on its mass and also the change of the orbital trajectories with 
changing orbital angular momentum. We use this simplistic ap- 
proach as it covers best the effect of a possible reaction of the 
host galaxy on the intruder (however, see Taylor & Babul ( 2001 ) 
and |Gan et aL] ( |2010[ ) for a more elaborated approach to model 
dynamical friction). 

Thus above a certain mass the satellite will not be able to eject 
any HVSs since the energy loss of the whole system is larger 
than the energy gain of the single stars. However, judging from 
our simulations this will only happen at masses > 10 11 M G . 
Extrapolating our results into this (major merger) mass regime 
is not meaningful as a massive intruder will significantly perturb 
the host galaxy. We can thus just state that this scenario does not 
occur for minor mergers in the present work. 

Hypervelocity stars In the context of hypervelocity stars which 
are unbound to the total system we have to compare AE with 
E sat , the orbital energy of the satellite, since the energy E of such 
an unbound object must be 



E = E sat - A£ DF + AE > 0. 



Thus the condition 



aAE m „ c > -E m + AE DF . 



(15) 



(16) 



must be fullfilled to allow unbound stars to be generated during 
such a satellite-host galaxy encounter with a 0.45. 

5. The energy distribution 

A more robust measure of range of orbital energies covered by 
the stripped-off stars can be obtained when looking on the over- 
all width of the energy distribution of all particles gravitationally 
unbound to the satellite. We determine these stars using an iter- 



ative method described by Tormen et al. (T998J). The energy dis- 
tribution is plotted in Fig. [6] The two peaks represent the leading 
and trailing tidal arm, respectively. The gap between them is cen- 
tered on the orbital energy of the satellite. We apply an empirical 
fitting function 



ME) 



/high 

c 



1 



(1 +exp(y(^ 



I))) 2 



./inner CXp 



AE 



2\ 



100 



80 



60 



40 



20 



— - Fitting function 




Trailing arm;l ;HVSs 



-40 -30 -E„ t -20 -10 
Orbital energy E [10 4 km 2 s~ 2 ] 

Fig. 6. A histogram of the orbital energies of all particles which 
initially belonged to the satellite but became unbound to it in the 
course of the first orbit (tidal debris particles). The two peaks 
correspond to the two tidal arms torn out of the satellite. The 
central gap coincides with the orbital energy E sat of the remain- 
ing satellite system. The dashed rad line shows the fitting func- 
tion (Eq. 17 1 which is used to characterize the distribution. The 
meaning of the fitting parameter e w which is used to determine 
the width of the high energy peak is also indicated. Note, that 
throughout this work we only consider the high energy peak, as 
this where HVSs would reside. 



again with AE = E - E sat . To normalize the function we use the 
factor /high/C, where /hi gn is the number of particles in the high 
energy (trailing) arm divided by the initial number of satellite 
particles and 



C = 



'Em 
1 + 



-(Hi 



7\ _ 



1 +exp( 



-7)) 



inner Mnner 



(18) 

The function is shown in Fig.|6]as red dashed line. For the fitting 
procedure we only consider the high energy peak of the distri- 
bution, i.e. where E > E s . dt . The fitting parameters provide us 
with some characteristics of respective distribution: the width or 
typical energy, e w , the width of the central minimum, ejnner. a 
measure of how fast the distribution drops off with increasing 
energy, y, and the relative depth of the central minimum, /inner- 



A composite distribution Following an idea by |Johnston] (l998) 
we assume the general shape of the distributions to be invariant 
and only the width and the normalization to be specific to the 
respective orbit and satellite parameters. This means that the pa- 
rameters ejnner, J, /inner are either constant for all situations or a 
function of e w . To compare the distribution shapes we rescaled 
the particle energies into units of their specific typical energy e w 
via 

Ei — -Esat 



AE-, 



(19) 



(17) 



The resulting energy distributions were then renormalized to 
eliminate the influence of the number particles in the respective 
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Fig. 7. Upper panel: The histograms of the energies differences Fig. 8. Relation between the maximum energy gain AE B 



AE = Ej-E s:i t of the tidal debris particles of all simulations used 
in this work overplotted. The distributions are plotted as func- 
tions of AE in units of the width of the high energy peak e w . The 
histograms are also renormalized so that the high energy peak 
covers the same area. The lines are color-coded according to the 
initial angular momentum of the progenitor system (from red 
representing more radial orbits to blue for more circular orbits) 
Lower panel: The mean distribution obtained from the distribu- 
tion plotted in the upper panel (solid black line). The standard 
deviation is also indicated with the thin grey lines. Overplotted 
in red is our fitting function using the parameters given in Eqs. 20 
and ew = 1 . 



(trailing) tidal arm. The upper panel of Fig. [TJoverplots the en- 
ergy distributions of all our simulations. Note, that these include 
satellite varying over more than a magnitude in mass and angu- 
lar momentum. The resulting mean distribution (black line) with 
the standard deviation (grey line) is also plotted. Applying our 
fitting function Eq. fTTlresults in the following parameter values: 



dinner — 0.146^ 

r = 5.21, 

Jhmer ~ 0.94. 



(20) 



The fit is shown as a dashed red line in Fig. [7] We then repeat 
the fitting procedure for all single distributions with only e w as a 
free parameter. The resulting widths are plotted in Fig. [8] against 
the highest energy of all satellite particles, AEnux.shn- Grey dots 
represent the lower resolution runs. The tight correlation 



AE„ 



■ft^ 



(21) 



with ft - 1.5 (which is not strongly affected by the resolution in 
the simulations) is a result of the steep drop at the high energy 
tip of the distribution. Using the simple model developed in the 
previous section we can thus obtain an estimate for the width of 
the high energy peak via 



£w — ft A.E'max.sim — ^ AE n 

— 0.3AE maxc . 



and the width e w of the energy distribution of all tidal debris 
stars obtain via fitting Equation 17 The grey dotted line has a 
slope of 1.5. 

Note, that v esc (r) is a proxy for the mass of the satellite. Thus 
more massive satellites will produce a larger energy spread in 
the stripped-off stars. One could also say that the higher velocity 
dispersion of a more massive galaxy directly translates into a 
larger energy dispersion in the tidal debris. 

6. Discussion 

It is now straight forward to compute the maximum velocities 
generated during a tidal collision at certain galactocentric radius 



<(r) = [2(£ max -<D h0St (r))]i 



[2(£ S at,a 



apo "t" 0-45AZ£max,c 



4> host (r))] 5 



(23) 



where E satApo — E sat - AE DF is the orbital energy of the remain- 
ing satellite after the passage. To obtain a quantitative idea we 
assume again a Galactocentric distance of r — 60 kpc to be com- 
parable to the observations of HVSs. For the energy loss via dy- 
namical friction A£df for simplicity we use the loss law (Eq. 14 1 
found in our simulations. The velocities obtained in this way are 
plotted in Fig.|9]as a function of the initial satellite mass M sat and 
for four different initial angular momenta L sat . The satellite sys- 
tem was assumed to approach the host galaxy on a parabolic or- 
bit, i.e. E sa t = km 2 s~ 2 . Only the most massive satellite galaxies 
eject HVSs with substantial velocities comparable to the radial 
velocities of the observed HVSs. Less massive galaxies could, in 
principal, also yield such large velocities if they move on more 
energetic orbits themselves. For example, a satellite with mass 
A^sat = 10 9 M would have to cross the virial radius of our par- 
ent galaxuHwith a velocity of ~ 660 km s _1 to eject a star with 
720 km s at 60 kpc comparable to the fastest HVSs known. 
Such a system would not lose enough energy to become bound 
to the larger galaxy. In this case one hardly would speak about 
an ejected star since the galaxy would move along with the star 



(22) 



In our Galaxy model ^200 - 260 kpc. 
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Fig. 9. The maximum ejection velocities at a galactocentric dis- 
tance of 60 kpc as a function of initial satellite mass as computed 
from Eq. 23 assuming an initial orbital energy of the satellite 
fisat = knrs~ 2 . The energy loss due to dynamical friction was 
computed using the empirical law (Eq. [T4] i obtained from our 
simulations. The lower grey lines show the velocity of the satel- 
lite remnant at the same distance. 



for a long period of time. 

We now consider the subsample of observed HVSs with travel 
times ~ 133 Myr pointed out by Abadi et al. (2009 cf. their 
Fig. 1). The spread in velocities is roughlj|j] 400 km s _1 . Such 
strong variations in the velocities translate into a progenitor mass 
of =i 10 11 M Q . If we select only stars within the overdensity re- 
gion defined by Abadi et al. (2009) the spread reduces to ~ 250 
km s _1 resulting in a minimum progenitor mass still larger than 
10 10 M . Concerning the known satellite galaxies the respective 
authors of the HVS discovery papers already excluded a kine- 
matical connection to them (|Brown et aL||2005| |2006a| [2009a; 
Edelman n et al.|2005| |Hirsch et al.|2005j ! Since we do not ex- 
pect such a system to have escaped observations to date we con- 
clude that a satellite origin for the subsample to be unlikely. 
However, one should keep in mind that a more massive host 
galaxy (our host model has a total mass of 1.1 x 10 12 M ) would 
shift the lines in Fig. [9] upwards and would, in principal, al- 
low also very small galaxies to produce stars with velocities, 
e.g. > 500 km s . However, the fact that only massive galax- 
ies can eject stars with velocities significantly larger than their 
own velocities remains unaffected by this. 



Fig. 10. Upper panel: Mass ejected at unbound velocities during 
one satellite orbit as a function of initial progenitor mass. The 
satellite is assumed to have zero initial orbital energy, i.e. to be 
on an parabolic orbit. Different lines correspond to four different 
initial angular momenta of the satellite. The vertical dashed line 
indicates the mass of a single star particle in our simulations. 
Lower panel: mass fraction of an unbound intergalactic popu- 
lation originating from satellites with masses below M sat . The 
fraction were computed using satellite mass function based on 
the satellite luminosity function of Koposov et al. (2008). The 
mass function is also shown as thin dashed line. More than 95 
percent of the population is created by satellites more massive 
than 10 9 Mq. 



This means that the phase space distribution of stellar halo stars 
reaches all velocities up to the local escape speed at all times. 
For example, |Smith et~aL ( 2007) 1 used this as a critical assump- 
tion for their technique to estimate the mass of the Galaxy. 
However, this also means that a classification of a star as a HVS 
ejected from the SMBH based on its velocity is only valid for ex- 
tremely large velocities. Without a confirmation of their young 
ages the "bound" HVS population in the compilation of [Brown 
|et al.| ( |2009a| l is indistinguishable from the normal (accreted) 
stellar halo population. To date only three stars in the survey 
have clear spectroscopic identification as main sequences B stars 



( Fuentes et al.|2006 Lopez-Morales & Bonanos 2008 Przybilla 



et al. 2008 1 while others could be old blue stragglers or blue 
horizontal branch stars ( |Perets et al.|2009| ). 



The bound HVS population and the outer stellar halo Several 
recent studies have shown that the outer stellar halo is almost 
purely made of accreted stars ( |Abadi et al.|2006| |Zolotov et al. 
2009 ; Scannapieco et al. 2009J. As smooth gas accretion via cold 
flows plays only a minor role for Milky Way-type galaxies at low 
red shifts (Brooks et al. 2009| l, we can assume the Galaxy has not 
grown significantly since its last major merger. Our simulations 
now demonstrate that satellite accretion will inevitably produce 
stars with velocities up to and exceeding their local escape speed. 



2 We ignore the fact that the stars reside at different galactocentric 
radii. Taking this into account would enlarge the spread even further. 



An intragalactic stellar population With our results we can also 
address the question on what kind of satellites are the main con- 
tributors to a possible intragalactic stellar population (ISP) or 
Wandering stars ( |Teyssier"et al. 2009). For this we assume the 
in-falling satellite galaxies to be initially on parabolic orbits, i.e. 
£sat = 0. Via dynamical friction the satellites will be shifted onto 
bound orbits during their first passage. Thus not all stars in the 
trailing (high energy) tidal arm will become unbound, but only 
those which gained more energy than is lost by their progenitor 
system. 



By integrating our fitting formula (Eq. 17 1 using the proper value 
for e w obtained via Equations [13] and |2 1 1 over energies larger 
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than the frictional energy loss A£df we obtain the fraction of 
the baryonic mass which became HVSs. In the upper panel of 
Fig. 10 this fraction multiplied by the baryonic mass content of 
the satellite is shown as a function of the total satellite mass. For 
the estimation an initial angular momentum Lq had to be set. The 
plot shows the results for four different Lq (color coding is the 
same as in Pig. |9j. 

We then convolve this mass ejection function with an obser- 
vational mass function of dwarf galaxies. We therefore used the 
luminosity function obtained by Koposov et al. (2008 ) and con- 
verted it into a mass function using the same relations used to 
create our satellite models for the simulations. As a result we 
obtained the cumulative HVS mass production function shown 
in the lower panel of Fig. [10] Color coding is the same as in 
the upper panel. The dashed black line indicates the cumulative 
mass function of satellite galaxies. Consistent with the results of 
|Teyssier et al. (2009) obtained via cosmological simulations we 
find that a tiny minority of massive satellites produces the over- 
whelming majority of HVSs. Given the fact that more massive 
galaxies usually are also more metal -rich (Tam ura et al.| [2001 ) 
we conclude that an intragalactic stellar population should have 
at least on average a higher metalicities than the surviving dwarf 
galaxy population. 

7. Summary 

In the present work we have used a suite of 41 A^-body simula- 
tions to study the tidal debris of satellite galaxies interacting with 
their much more massive host systems. Abadi et al. (2009) sug- 
gested that a fraction of the stripped-off stars can reach signifi- 
cant velocities and could be confound with Hypervelocity stars 
(HVSs) ejected from the Galactic center by a supermassive black 
hole. We find that, as suggested by these authors, the stripped- 
off stars are in fact observed in a confined region on the sky. 
However, for stars at distances still observable from the solar 
position the reason for this is not only the projection of a colli- 
mated stellar stream along the line of sight, but in addition that 
so shortly after the stripping event the stars had not yet enough 
time to disperse in physical space. 

We further developed a simple analytic model to predict the 
maximum possible ejection velocities via estimating the maxi- 
mum possible energy kick a star can obtain during such a tidal 
encounter (Eq. [13). 



Convolving our formalism with a satellite mass function al- 
lows us to determine the masses of the progenitors of the main 
contributors to a potential intergalactic stellar population (ISP). 
We find that stars originating from satellite galaxies with masses 
> 10 9 M Q form about 95 percent of the population. This is con- 



Following Johnston ( 1998 1 we suggest that the general shape 
of the energy distribution of particles stripped-off during one 
orbit is self-similar and can be described quite accurately by 
Equation 17 There are only two free parameters in the distribu- 



tion, its width and its normalization. The first represent a charac- 
teristic energy and is tightly connected to the maximum energy 
kick described by our stripping model. The normalization sim- 
ply reflects the fraction of mass lost by the satellite. Both can 
be predicted knowing only the initial properties of the host and 
satellite galaxy without the need of computationally expensive 
A^-body simulations. 

We also address the recently reported Hypervelocity star 
population. Velocities larger than 500 km s _1 are only gener- 
ated by massive satellite galaxies (> 10 11 M ) or by galaxies 
with very large infall velocities in which case these galaxies 
stay unbound from the host and leave the parent galaxy together 
with the HVSs. Furthermore the larger spread in velocities of 
HVSs with common traveltimes also requires a massive progen- 
itor (> 10 10 M Q ). The absence of the remnant of such a massive 
system makes a tidal debris origin for the HVSs unlikely even 
from a kinematical point of view. 



sistent with the findings of Teyssier et al. ( 2009 1 who traced back 



the origin of unbound particles in the cosmological simulations 
of Bulloc k~& Johnston] ([2005). We conclude thus that such an 
ISP should tend to have the same or even a higher metallicity 
than the outer halo population and also as the present population 
of Milky Way satellites. 
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Appendix A: Estimating the maximum energy gain 

In this section we briefly describe the process of how we ob- 
tained estimates for the velocities y per i and v esc (/? t id a i) which we 



need to evaluate Equation 12 The ingredients for this are 



- the radial mass profile of the host galaxy, 

- the radial dark matter and baryonic mass profile of the satel- 
lite galaxy, 

- the parameters of the satellite orbit, namely the initial angu- 
lar momentum L sat and the initial orbital energy E sat . 

In a first step we estimate the minimum distance to which the 
satellite approaches the host center, i.e. the pericenter distance 
^perf- For this we use the effective potential <t> e ff = ^g at / (2r 2 ) + 
<D host (r) exploiting 

£sat - ^AEdf = O e ff(7?p e ri), 



(A.l) 



where we compute the energy loss from dynamical friction AEdf 
using Eq. [14] Further we compute the satellite velocity in the 
Perigalacticon via 



V, 



peri 



2(£ sat - ^A£ DF - O host (tf peri )). (A.2) 



To compute the escape velocity v esc (^tidai) from the satellite 
system we first have to determine the tidal radius ,R t idai which we 
assume to be equal to the Jacobi radius at the distance R per i'. 

* tldd = l3M host (tf perl) ) (A3) 

However, we do not take the total satellite mass M sat for the final 
radius. We also take into account that due to its much larger ex- 
tension the dark matter halo of the satellite is stripped much ear- 
lier the the baryonic component. Consequently in a first step we 
compute the tidal radius using the total satellite mass M sat and 
assume that all material outside this "dark matter tidal radius" 
^udai.DM is lost. We then compute the "baryonic tidal radius" us- 
ing Equation A. 3 with the mass M' sat - M sat (r < ^ t id a i,DM)- 



Finally we obtain the escape speed 

1 esc (#tidal) = V2|<t> sat (fltid a l)l 



(A.4) 



Note, that the tidal radius computed in this two-step process also 
allows a very good estimate of the baryonic mass loss of the 
satellite when it is assumed that all mass outside this tidal radius 
is lost, i.e. 



Most M sat (r > fl tidal ) 



unbound 



(A.5) 



M sat M sat 

This was used in Sect.|6]to estimate the fraction of satellite mass 
expelled as HVSs into intergalactic space. 
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Fig. B.l. Comparison of the energy distributions obtained from 
corresponding high and low resolution runs. The dashed lines 
indicate the mass resolution limits of the simulations, i.e. the 
mass of a single star particle. 



Appendix B: Scaling tests 

To assess the influence of the numerical resolution on our re- 
sults we ran a set of simulations with the number of particles 
in the satellite system reduced by a factor of 10. This was done 
by taking the initial snapshot file of one of our high resolution 
runs and randomly removing 90 percent of the satellite particles. 
Then the masses of the remaining particles increased by a factor 
of 10. In this way we obtain an equilibrium configuration of a 
satellite system with exactly the same properties as the high res- 
olution one. 



In Figure B.l we show the resulting energy distributions to- 
gether with those of the corresponding high resolution runs with 
the same initial conditions. Dashed lines indicate the respective 
mass resolution limits, i.e. the mass of a single star particle in 
the simulations. The shapes of the distributions show no signif- 
icant differences. The width of the corresponding distributions, 
e w , obtained by fitting Eq. 17 also differ by less than 10 percent. 
As could be expected the low resolution distribution does not 
reach as high energies as the high resolution one. However, the 
maximum energies do not differ by much, due to the steep slope 
in the outer tails of the distribution. 

We also repeated one simulation run using a 5 times larger 
softening length for the star particles. For all quantities mea- 
sured for this study the outcome changed by less than 1 percent. 
Especially the maximum energies reached by satellite particles 
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Fig. C.l. Radial gravitational potential profiles of the four alter- 
native host galaxy representation and for the A^-body live host 
galaxy (black). The alternative models cover a variety of central 
and outer slopes allowing to test the influence of those on satel- 
lite tidal tails. The varying thickness of the profile lines reflects 
the non-spherical symmetry of the potentials. 
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rigid potentials representing the host galaxy. Four different mod- 
els were applied: three of them (potentials 1-3) share the same 
pa rametrization of the barvon ic disk and spheriod components. 
A |Miyamoto & Nagai] ( |l975| l disk with mass 10 11 M , radial 
and vertical scale length of 6.5 and 0.26 kpc, respectively and a 
Hernquist ( 1990) bulge component of mass 3.4 x 10 10 M and a 



scale length of 0.7 kpc. The dark haloes are modeled with 

- potential 1: a logarithmic potential O(r) = v? ln(r 2 + r^) with 
v'h = 128 km/s and rj, = 12 kpc. This potential was already 
often used by other studies to mimick a Milky Way (e.g. 
|Johnston| 19981 IHelmi & White|2001) . 

- potential 2: a |Plummer| ( |191 1| sphere with mass 6 X 10 11 M G 
and scale radius 25.7 kpc. 

- potential 3: a NFW sphere ( |Navarro et al. 1997| l with central 



density 1.523 x 10 6 M kpc and scale radius 36 kpc 

All three potentials were chosen such that V c irc(8.5 kpc) = 220 
km/s and V esc (8.5 kpc) = 550 km/s. As a fourth option we 
used the model potential number 4 of |Dehnen & Binney ( 1998) 
which we implemented into Gadget-2 using a C++ routine pre- 
pared by Walter Dehnen and distributed with the NEMO Stellar 
Dynamics Toolbox (Teub en|1995 1. 

Figure C.l shows a comparison of the radial profiles of the 
four potentials with the radial profile of the live halo used in the 
main part of the simulations. All potentials have a steeper slope 
in the inner regions. The virtually flat part of the live potential is 
due to the gravitational softening becoming significant on these 
scales. 

Figure [C2l plots the energy tidal tail distributions (cf. Fig. [6j ob- 
tained with the different host representations but otherwise iden- 
tical initial conditions. While the distribution changes strongly 
in regions with small \AE\, the tail of the distribution remain 
virtually unchanged. We thus conclude that the actual shape of 
the Galactic potential has no major influence on our results. The 
variations around the central minimum are most likely due to 
different evolution of the Roche radius of the satellite during its 
orbits determining whether particles with low |AE| which stay 
near the satellite for longer periods are re-captured. 



Appendix D: Initial conditions 



Fig. C.l. Comparison of the energy distributions of the tidal tail 
particle stripped from identical satellite galaxies with identical 
initial phase space positions evolving in different host potentials. 
The outer tails are virtually unaffected by the differing potentials, 
while the central minimum is subject to significant changes. 



differ only by 0.1 percent. This shows that our results are not 
affected by artificial heating by two-body encounters. 



Appendix C: Other host galaxy potentials 

To test the influence of the shape of the Milky Way-host potential 
we performed a small number of test simulations with different 
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Tilmann Piffl et al.: Kinematics in galactic tidal tails 
Table D.l. Initial parameters of the satellite systems (plus some analysis results) 



No. 


N 1 , 


Km 


200 


M 4 

barv 




c o 


orbit type 7 


4 


Jul) 










10 10 M 


10 8 M o 


10 3 kpc km s~' 


10 4 km 2 s~ 2 




10 4 km 2 s~ 2 




Runs with reduced resolution 


001 


10 4 


5 x 


10 3 


1.0 


1.0 


0.5 


-3.8 


polar 


2.2 


0.49 


002 


10 4 


5 x 


10 3 


1.0 


1.0 


3 


-1.8 


polar 


1.9 


0.30 


003 


10 4 


5x 


10 3 


1.0 


1.0 


4 


-0.9 


polar 


1.8 


0.15 


004 


10 4 


5 x 


10 3 


1.0 


1.0 


7.5 


-1.8 


polar 


1.6 


0.10 


005 


10 4 


5 x 


10 3 


1.0 


1.0 


15 


-1.5 


polar 


1.7 


0.05 














Rune; mi tVi ■fill 1 


tp <^n1 lit i nn 








006 


i a5 


5 x 


a4 

10 


0.7 


0.8 





-1.8 


polar 


2. 1 


0.65 


007 


10- 


5 x 


a4 

4 


1.0 


1.0 





-1.8 


polar 


2.6 


0.56 


008 


10 


5 x 


a4 

10 


0.1 


0.06 


0.5 


-3.8 


polar 


0.6 


0.71 


009 


10- 


5 x 


a4 

10 


0.3 


0.2 


0.5 


-3.8 


polar 


1 . 1 


0.58 


010 


i 

10 


5 x 


a4 

10 


0.5 


0.5 


0.5 


-3.8 


polar 


1.6 


0.55 


011 


i a^ 

10- 


5 x 


10 


0.5 


0.5 


0.5 


-3.8 


polar 


1.6 


0.55 


012 


i AS 
10 


5 x 


a4 
4 


0.7 


0.8 


0.5 


-3.8 


polar 


1.9 


0.51 


013 


i a^ 

10- 


5 x 


a4 

K) 4 


1.0 


1.0 


0.5 


-3.8 


polar 


2.4 


0.47 


014 


i AS 
10 


5 x 


a4 
4 


2.0 


3.0 


0.5 


-3.8 


polar 


3.7 


0.42 


015 


i a^ 

10- 


5 x 


rA 




0. 1 


0.06 


3 


-1.8 


polar 


0.7 


0.26 


016 


i 

10 


5 x 


a4 

10 


0.3 


0.2 


3 


-1.8 


polar 


1.1 


0.23 


017 


i a^ 

10- 


5 x 


a4 

K) 4 


0.5 


0.5 


3 


-1.8 


polar 


1.5 


0.21 


018 


10 5 


5 x 


a4 

10 


0.7 


0.8 


3 


-1.8 


polar 


1.7 


0.20 


019 


i a^ 

io- 


5 x 


a4 


1.0 


1 .0 


3 


-1.8 


polar 


2.0 


0.19 


020 


i 

10 


5 x 


a4 

10 


1.0 


1.0 


3 


-1.8 


polar 


2.0 


0.20 


021 


i a^ 

10- 


5 x 


a4 

K) 4 


2.0 


3.0 


3 


-1.8 


polar 


2.9 


0.20 


022 


10 5 


5 x 


a4 

10 


1.0 


1.0 


3 


-1.8 


polar 


2.0 


0.19 


023 


i a^ 

io- 


5 x 


a4 

10 


0.3 


0.2 


15 


-1.5 


polar 


1 .0 


A CI/. 

0.06 


024 


i 

10 


5 x 


a4 

10 


0.5 


0.5 


15 


-1.5 


polar 


1.2 


0.06 


025 


i a^ 

io- 


5 x 


a4 
.0 


0.7 


0.8 


15 


-1.5 


polar 


1 .4 


A A£ 

0.06 


026 


i 

10 


5 x 


a4 

10 


1.0 


1.0 


15 


-1.5 


polar 


1.6 


0.05 


027 


i a^ 

io- 


5 x 


a4 

10 


0.5 


0.5 


3 


-2.8 


polar 


1 A 

1 .4 


a n 

0.23 


028 


10 


5 x 


A4 


0.5 


0.5 


3 


-1.2 


polar 


1.4 


a in 

0.20 


029 


i a^ 

10- 


5 x 


a4 

K) 4 


0.5 


0.5 


3 


-0.8 


polar 


1 .4 


0.20 


030 


i aS 
10 


5 x 


a4 

10 4 


1.0 


1.0 


7.5 


-1.7 


polar 


1.8 


0.10 


031 


i a^ 

10 


5 x 


rA 
10 


1 .0 


1.0 


7.5 


-1.7 


45° retrograde 


1.8 


0.10 


032 


i aS 
10 


5 x 


a4 

10 4 


1.0 


1.0 


7.5 


-1.7 


45° prograde 


1.8 


0.11 


033 


1 TO 

10 


5 x 


rA 

10 


1.0 


1.0 


3 


-1.8 


planar retrograde 


2.1 


0.25 


034 


i aS 
10 


5 x 


a4 

10 4 


1.0 


1.0 


3 


-1.8 


planar prograde 


2.0 


0.25 


035 


10 5 


5 x 


10 4 


1.0 


1.0 


7.5 


-1.7 


planar retrograde 


1.8 


0.11 


036 


10 5 


5x 


10 4 


1.0 


1.0 


7.5 


-1.7 


planar prograde 


1.8 


0.12 


037 


10 5 


5x 


10 4 


0.5 


0.5 


15 


-1.5 


planar retrograde 


1.2 


0.06 


038 


10 5 


5x 


10 4 


0.5 


0.5 


15 


-1.5 


planar prograde 


1.2 


0.06 


039 


10 5 


5x 


10 4 


1.0 


1.0 


1 


-2.8 


polar 


2.6 


0.33 


040 


10 5 


5x 


10 4 


1.0 


1.0 


1.6 


-2.4 


polar 


2.3 


0.27 


041 


10 5 


5x 


10 4 


1.0 


1.0 


2 


-2.4 


polar 


2.2 


0.25 


042 


10 5 


5x 


10 4 


1.0 


1.0 


4 


-0.9 


polar 


1.9 


0.15 


043 


10 5 


5x 


10 4 


0.5 


0.5 


7.5 


-1.8 


planar retrograde 


1.4 


0.13 


044 


10 5 


5x 


10 4 


1.0 


1.0 


1 


0.0 


82° prograde 


2.5 


0.42 


045 


10 5 


5x 


10 4 


1.0 


1.0 


5 


0.0 


87° retrograde 


2.0 


0.12 


046 


10 5 


5x 


10 4 


1.0 


1.0 


10 


0.0 


84° retrograde 


1.9 


0.06 



1 number of star particles in the satellite 
2 number of particles in the satellite dark halo 
3 virial mass of the satellite 
4 mass of the satellite baryonic component 
5 initial angular momentum of the satellite with respect to the host galaxy 
6 initial orbital energy of the satellite 
'inclination of the satellite orbital plane with respect to the host disk plane and, 
if appropriate, sense of rotation relative to the host disk rotation 
8 width of the resulting energy distribution as defined in Eq. 
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fraction of baryonic mass lost during the first orbit of the satellite. 



